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Abstract 

The  objectives  of  this  stiiUy  are  to  characterize  apparent 
background  motion  as  a  function  of  observer-satellite  orbit  and 
sighting  geometry  aiul  then  to  compare  the  effectiveness  of  signal 
processing  filters  in  minimizing  background  motion  effects  on 
target  detection.  I'or  many  applications,  complex  algorithms 
cannot  he  realized  economically,  so  simple  yet  effective  algorithms 
must  be  developed  to  do  the  job.  In  this  study  both  ideal  and 
optimum  third-order  temporal  recursive  filters  were  synthesized 
and  their  simulated  signal-to-noise  ratio  perfomtance  compared 
with  that  obtained  by  simple  temporal  differencing  (first-  through 
fourth-orderl.  Our  results  indicate  that  for  the  same  amount  of 
memory,  compared  to  the  differencing  filters,  significant  SNR 
performance  improvement  can  be  achieved  with  optimum  recursive 
filters  at  the  cost  of  a  modest  increase  in  filter  complexity.  The 
study  is  particularly  pertinent  to  the  detection  by  a  staring  mosaic 
sensor  of  aircraft  from  satellites  in  high  (>  1000  km  I  orbits.  They 
are  typically  dim  targets  fiying  close  to  a  structured  background 
(the  earth  surface I.  and  the  background  clutter  tends  to  dominate 
the  detection  problem. 

I.  Introduction 

Staring  mosaic  sensors  are  effective  in  minimizing  possible 
"noise"  from  background  spatial  structure  by  taking  differences 
between  successive  frames  of  data  and  thereby  eliminating  features 
which  have  spatial  but  not  temporal  variation.  This  is  commonly 
referred  to  as  background  subtraction.  Thus  in  space  surveillance, 
the  star  background  can  be  eliminated  by  an  inertially  stabilized 
staring  sensor  while  "targets"  (satellites  which  move  through  the 
field  of  view  I  are  detected  as  streaks.  Detection  of  objects  against 
an  earth  background  is  more  difficult  because  the  staring  feature 
cannot  be  ideally  achieved.  In  any  orbit  except  synchronous 
equatorial,  an  earth-looking  sensor  can  effect  perfect  staring  at 
only  one  point  in  the  FOV  because  of  both  satellite  motion  and 
earth  motion.  All  other  background  points  in  the  FOV  will  appear 
to  drift,  and  this  motion  causes  "noise"  signals  if  the  background 
has  appreciable  spatial  structure  to  simulate  one  or  more  tem¬ 
porally  moving  targets.  Even  in  a  synchronous  equatorial  orbit, 
attitude  drift  of  the  sensor  is  of  concern  in  establishing  ideal  staring 
operation.  The  "expansion"  and  "contraction"  of  the  FOV, 
depending  on  the  satellite  zenith  angle,  is  shown  in  Fig.  I  (earth 
curvature  over  the  satellite  FOV  is  ignored). 

A  critical  problem  in  the  design  of  signal  processing  systems 
for  large  staring  mosaic  (detector  array)  sensors  is  the  synthesis  of 
algorithms  to  detect  dim  targets  against  earth  backgrounds  in  the 
presence  of  several  noise  sources.  Significant  noise  components  arc 
background  drift,  sensor  jitter,  background  shot  noise,  detector 
crosstalk,  system  readout,  and  multiplexer  (MUX)  noises.  In  many 
applications  the  implementable  signal  processing  algorithms  are 
limited  by  the  large  number  of  detector  channels  required  for 
adequate  spatial  sampling  resolution  and  coverage.  Consequently, 
filtering  algorithms  must  be  simple  and  must  provide,  under  given 
constraints  on  memory  and  computing  capacity,  reliable  detection 
at  low  input  signal-to-noise  ratios  against  severe  and  varying 
background  conditions.  The  overall  performance  of  a  dim  target 
detection  system  is  limited  by  the  amount  of  SNR  enhancement 
I  which  can  be  achieved  by  a  filter  applied  to  the  raw  focal  plane 

data.  Because  of  the  large  amount  of  data  generated  by  a  mosaic 
sensor,  a  primary  goal  of  each  detection  algorithm  block  is  to 
reduce  the  dimensionality  of  the  problem  such  that  the  final  step 
*  -■*  can  be  a  simple  threshold-type  classification.  To  first  order,  sensor 

^  ^  transfer  characteristics  are  linear  as  shown  in  Fig.  2.  Ideal  and 

optimum  filters  can  be  synthesized  directly  from  linear  models  of 
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the  sensor  and  its  inputs.  In  the  context  of  this  report,  an  "ideal" 
filter  is  one  that  is  not  constrained  to  be  realizable  while  an 
“optimum"  filter  must  be  tealizable  within  the  constraints  of  a 
given  filter  configuration.  This  paper  will  include. 


Discussion  of  the  expected  earth  background  drift 
velocities  that  a  staring  sensor  in  a  low  earth  orbit  (LEO) 
will  encounter,  and  characterization  of  such  drift  in 
terms  of  orbital  and  earth  viewing  parameters. 
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I*rosonl;iiion  ol' simpit.-  temporal  filtering  algorithm!;,  anti 
a  comparative  analysis  of  the  effectiveness  of  various 
signal  processing  algorithms  in  rejecting  background 
clutter  noise. 

I  he  interrelationship  between  blur  circle  size  and  the 
consetpient  theoretical  limit  in  achievable  SNR. 

llic  results  summari/cd  in  this  paper  were  obtained  under  a 
Rockwell  IR&l)  effort  and  are  e.xplained  in  more  detail  in 
Reterence  1.  Some  earlier  results  obtained  under  this  same  IR&D 
etiort  were  reported  in  Reference  2  and  also  are  included  in  this 
paper. 


the  motion  of  any  point  in  the  sensor  FOV  is  linearly  proportional 
to  its  distance  from  the  FOV  center,  the  stable  point  on  the  ground. 
Thus,  results  shown  in  this  study  for  the  perimeter  of  a  square 
FOV  can  be  extrapolated  to  larger  FOV’s  or  to  circular  FOV's. 

System  parameters  used  to  study  background  motion  depend¬ 
ence  on  sighting  geometry  are  shown  in  Table  I.  For  the  first  part 
of  the  study  to  determine  the  dependence  of  drift  parameters  on 
satellite  orbit  altitude  and  zenith  angle,  three  orbit  altitudes- 
400n.mi.  (741.2km).  lOOOn.mi.  (1853  km).  and  3440n.mi. 
(6374  km)-were  picked  and  for  each  orbital  altitude,  four  in-plane 
sighting  geometries  were  chosen  with  satellite  zenith  angles  of  30*’. 
45°.  60°  and  80°. 


II.  Characterization  of  Farth  Background  Drift 

It  we  ignore  the  earth's  rotation,  the  apparent  angular  motion 
ot  a  ground  point  on  the  edge  of  the  sensor  FOV  is  given  by  ^ 


Vs 

cJ  = -  xAro,. 

r|  ro  ‘Js 


where,  referring  to  Fig.  3a. 

is  the  sensor  linear  velocity  in  the  satellite  orbit  plane 

7|  is  the  vector  to  the  point  of  interest  on  the  edge  of  FOV 

Tq  is  the  vector  to  the  F'OV  stable  point  on  the  ground 

Tqc  ^  vector  equal  in  magnitude  to  Ar 

A  similar  derivation  is  applicable  to  the  study  of  apparent 
motion  of  other  points  on  the  edge  of  the  sensor  FOV.  The  net 
motion  along  the  xyz  axes  with  the  x  axis  parallel  to  the  satellite 
velocity  in  the  orbital  plane  can  be  computed  by  resolving  the7| 
vector  into  three  components  (viz..?)^  Tjy  andTjy..)  and  then 
taking  tlie  vector  cross-product  of  eacli  with  the  satellite  velocity 
vector  Vj.  Fig.  3b  illustrates  such  a  procedure.  Depending  on  the 
position  of  the  point  of  interest  (target)  in  the  sensor  FOV,  it  can 
have  simple  linear  motion  or  combined  linear  and  rotational 
motion  about  the  FOV  center.  For  small  FOV’s,  the  magnitude  of 


Table  I.  System  Parameters  to  Study  Background  Motion 
Dependence  on  Sighting  Geometry 


Psrimtttr 

Cliaractarittia 

Circular  catallita  orbit 

Orbit  inclination 

0°  (equatorial) 

Orbit  altitude 

254  km  (137  n.mi.)  ■  6374  km  (3440  n.mi.l 

Satellite  orbit  period  range 

1.495  hr -4.00  hr 

Satellite  zenith  angle  WAT  local 

30-80® 

vertical  of  earth  aimpoint 
£arth  aimpoint/cet>teT  of  FOV 

Latitude  range  »  0-19^ 

Longitude  range  =  3.6-50.5° 

Longitude  at  apogee 

0“ 

Argument  of  perigee 

1/2  focal  plane  angle 

0“ 

Forward  FI 

Side  SI 

Sighting  geometry 

Along  orbital  plane  forward  of  direction  of 
motion:  60°  to  tide  in  plane  perpendicular 

to  orbit  plane 

There  are  two  more  meaningful  measures  of  angular  drift 
other  than  the  maximum  angular  drift.  lTO|  the  latter 
occurring  at  one  or  two  of  the  four  comers  of  a  square  or 
rectangular  FOV.  Using  linear  interpolation  it  is  easy  to  derive  the 
average  drift  (average  of  eight  points)  on  the  circumference  of  a 
circular  FOV  with  radius  r^,  jTO|  re.  circumscribed  inside  the 
square  FOV  as  in  Fig.  3c.  Given  the  value  of  |TO|  re-  h  is  simple  to 
derive'  the  value  of  the  second  important  drift  averaged  over  the 
area  of  a  circular  FOV  of  radius  r^  to  be 

lT0|F0V  =  f  |TO|re 


The  values  of  maximum  angular  and  linear  background  drift 
parameters  have  been  plotted  against  various  pertinent  orbital 
parameters.  Fig,  4  shows  the  magnitude  of  the  maximum  total 
angular  drift  iTOlmax  (arc-sec/sec)  as  a  function  of  the  satellite 
orbit  altitude  (Fig.  4a)  and  zenith  angle  (Fig.  4b).  The  maximum 
magnitude  of  the  angular  drift  usually  occurs  at  one  or  two  of  the 
four  comers  of  a  rectangular  FOV  and  has.  for  constant  zenith 
angles,  a  somewhat  inverse  relationship  to  orbit  altitude.  For 
constant  satellite  altitude,  the  magnitude  of  the  maximum  angular 
drift  lTO|  max  increases  markedly  with  higher  satellite  zenith 
angle,  although  the  steepest  increase  takes  place  for  satellite  zenith 
angles  greater  than  60°.  as  shown  in  Fig.  4b. 


Fig.  5  shows  the  various  meaningful  measures  of  angular  drift 
and  their  dependence  on  satellite  orbit  altitude.  Two  sets  of  curves 
are  plotted  for  the  same  satellite  zenith  angle  (viz.,  60°).  One  set 
lefers  to  line  of  sight  (LOS)  60°  forward  looking  in  the  orbital 
plane  and  the  other  set  refers  to  the  LOS  looking  60°  to  the  side  at 
right  angles  to  the  orbital  plane.  Typically,  the  angular  drift  rate 
looking  .sideways  is  significantly  less  than  looking  along  the  orbital 
plane  for  small  (<1000  km)  satellite  orbit  altitudes.  Average  drift 
magnitude  on  the  circumference  of  a  circular  FOV  is  roughly  70'; 
of  the  maximum  drift  value.  It  has  been  analytically  shown  that 
the  angular  drift  averaged  over  the  area  of  a  circular  F'OV  is  2  3 
that  of  the  angular  drift  on  the  circumference  of  the  same  FOV. 
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Fig.  4.  Dependence  of  Maximum  Total  A  ngular  Drift  on  Satellite 
Altitude  and  Zenith  Angle 
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The  variation  of  the  apparent  background  linear  drift  velocity 
orthogonal  to  LOS  with  orbital  altitude  is  less  than  one  might 
expect  because  angular  drift  has  an  inverse  relationship  with  slant 
range  for  the  same  satellite  zenith  angle.  This  is  illustrated  by  the 
two  curves  of  Fig.  6.  which  sliow  linear  drift  to  be  much  less 
dependent  on  orbital  altitude  over  the  range  shown.  This  is 
particularly  true  for  the  forward-viewing  case,  which  is  also  the 
worst  case. 


Fig.  5.  Dependence  of  Apparent  Backpound  Angular  Drift  Rater  on 
Satellite  Altitude  (Zenith  Angle  Mfi) 
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Flg.6.  Dependence  of  Apptrent  Background  Linear  Dr0 
Velocity  on  Satellite  Altitude 


III.  IVrl'oniiaiK'c  of  Various  Digital  Filters 
in  Suppressing  Background  Clutter 

One  approach  to  handling  background  motion  is  to  estimate 
this  motion  and  then  compensate  the  received  image  data  so  as  to 
remoNc  it.  The  algorithms  required  to  perform  high  accuracy  image 
motion  compensation  are  not  simple  enough  to  be  practically 
.ipplicabic  to  large  array  staring  mosaic  sensors.  Another  approach 
prewnted  here,  is  to  use  SNR  enhancement  filters  that  exhibit 
Millicicnt  perlormance  in  spite  of  excessive  background  motion. 
t»i\en  the  target  signal  magnitude  of  Fig.  7  and  sampled  back¬ 
ground  clutter  of  Fig.  8.  it  is  possible  to  use  a  high  band  pass  filter 
to  reject  the  low-frequency  background  clutter  and  retain  the 
high-frequency  components  of  signal  magnitude,  thus  significantly 
improving  the  input  SNR  throughput.  For  detecting  point  source 
targets  moving  against  structured  background,  we  use  primarily  the 
clas.sical  detection  ti.e..  pre-whitened  matched)  filttr.  This  filter 
can  he  written  as  H|(Ja))  =  S(-jtJ)/ |N(joJ)|  -  where  Sfjcu)  is  the 
signal  and  Ntjw)  the  additive  noise.  This  filter  optimizes  the  peak 
signal  to  rms  noise  ratio  (SNR)  and  is  typically  the  best  when  SNR 
gain  is  the  overriding  consideration.  In  Reference  2  it  has  been 
shown  that,  as  a  least-square  estimator,  the  realizable  filter  H^^fjea), 
as  an  approximation  to  the  ideal  filter  H|(juj),  is  optimum  with  an 
acceptable  delay  a  if  the  performance  parameter  given  by 

P=  f^lRljcjlj-  |H|(itj)-e"J“  HA(jcj||  -doj  (3) 

has  been  minimized,  where  R(Jcj)  is  the  Fourier  transform  of  the 
total  filter  input.  The  examples  of  the  optimum  Jinear  filter 
approximations  presented  in  later  sections  were  obtained  by 
minimizing  this  performance  parameter. 


Fig.  7.  Sampled  Signal  Magnitude  Spectrum  (l-Second  Sampling  Interval) 


Fig.  S.  Sampled  Noise  Magnitude  Spectra  (l-Second  Sampling  Interval) 


Because  of  the  nature  of  the  multiplexed  data  stream  out  of  a 
mosaic  sensor,  the  required  linear  filters  are  sampled-data  rather 
than  continuous.  Hence,  account  must  be  taken  of  fold-over  or 
aliasing  and  the  bandwidth  limit  specified  by  the  sampling 
theorem. 

In  the  simplest  approach  to  temporal  filtering,  one  matches 
the  target  detector  cell  dwell  time  to  the  detector  integration  time 
and  the  sampling  period  and  then  performs  some  order  of  temporal 
differencing.  A  slightly  more  complicated  approach  is  to  employ  a 
temporal  recunive  filter  and  sample  with  a  period  less  than  or 


equal  to  the  target  cell  dwell  time.  As  will  be  shown  later,  this 
second  approach  offers  greater  system  flexibility  and  performance 
at  a  slight  increase  in  complexity. 

To  compare  the  performance  of  various  filters  as  applied  to 
the  suppression  of  background  clutter,  a  computer  program  was 
used  to  simulate  a  set  of  fixed  (first  through  fourth)  difference 
filters  as  well  as  to  optimize  a  set  of  third-order  recursive  filters  for 
six  different  drift  rates.  To  show  that  the  fixed  difference  filters 
are  special  cases  of  optimum  filters  of  the  same  order,  we  start  with 
a  temporal  filter  with  a  third-order  recursive  structure  as  shown  in 
Fig.  9.  The  filter  is  assumed  to  have  a  transfer  function  H(Z) 
expressed  in  the  Z  domain  as. 


H(Z)  = 


Z  +  Oq 
Z  +  “p 


X 


Z-i-Kv 

Z-tK) 


Z2  +  2  coNo  Z  +  ‘^N0' 
Z2  +  2  JpCJNp  Z  +  tONp- 
Z2  -t  K4  Z  +  Kft 

z2  -h  K3  z  -t  K5 


(4) 


By  setting  O:o~l^2"'''  £  0  “  K4  =  -2  an.itojqO*” 
K5  =  I,  and  all  other  coefficients  equal  to  zero.  H(Z)  of  Fig.  9 
becomes  a  third  differencing  filter  H(Z)  =  (Z-l  )^/Z^)  which  is  a 
concatenation  of  three  first  difference  filters  [vis.,  (Z-1  l/Z] . 
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Fig.  9.  General  Third-Order  Recursive  Temporal  Filters 

Letting  all  of  the  filter  coefficients,  o  q-  iO-“NO-“D'  t 

and  WNp  of  Equation  4  vary,  a  set  of  six  optimum  thirdSirder 
recuisive  filters  was  synthesized  to  best  match  the  ideal  filter 
transfer  functions  according  to  the  criterion  set  by  Equation  3.  The 
coefficients  of  these  optimum  third-order  recuisive  filters  which 
depend  both  on  background  drift  rate  and  sampling  interval  are 
shown  in  Tables  2  and  3. 


# 
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Table  2.  Coefficients  for  Optimum  Recursive  Filters-Sampling 
at  Intervals  of  1/2  the  Target  Cell  Dwell  Time 


Vb 

Ipi/sl 

“0 

^0 

^MO 

ap 

“Ng 

7.3 

1.086 

-1.040 

0.756 

0.005 

0.091 

0.823 

6.1 

1.716 

-1.163 

1.758 

0.279 

0.146 

0.842 

4.8 

1.001 

6.612 

2.623 

0.253 

0.077 

0.824 

3.7 

-1.002 

2.323 

-1.379 

-0.193 

-0.367 

0.759 

2.4 

-0086 

1.000 

0094 

0.789 

0.678 

-0.774 

1.2 

0.205 

1.000 

1.015 

-0.370 

0.450 

-0.606 

Table  J  Coefficients  for  Optimum  Recursive  Filters  Sampling 
ai  the  Target  Cell  Dwell  Time 


Vh 

IptM 

"0 

^0 

“MO 

^P 

“Np 

7.3 

0.996 

-1.284 

0.287 

0.802 

0.768 

0.876 

6.1 

-1.000 

-0.854 

1.127 

0.680 

0.656 

0.714 

4.8 

-0.999 

-0.942 

0.884 

0.414 

0.464 

0.648 

3.7 

-1.003 

-1.004 

1.23 

0.462 

0.156 

0.704 

2.4 

-0.711 

-1.000 

0.963 

0.186 

0.448 

-0.706 

1.2 

-1.027 

-1.002 

1.055 

0.436 

0.920 

-0.411 
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resulting  spectra.  Fig.  8  is  especially  useful  as  it  also  illustrates  the 
effect  of  background  motion  on  background  clutter  spectra  and 
hence,  the  total  noise  spectra.  Fig.  8  is  indeed  a  super-position  of 
six  cases  for  six  background  drift  rates  (7.3p  r/s  to  1.2/ur/s  in 
equal  decrement  interval). 

Fig.  10  is  the  impulse  response  of  a  first  difference  filter.  The 
time  of  the  impulse  has  been  arbitrarily  set  at  SO  sec.  No'e  that  due 
to  the  property  of  first  difference,  we  obtain  two  samples  from  an 
unblurred  point  target,  of  equal  but  opposite  magnitude  and  set 
apart  by  the  sampling  interval.  The  impulse  responses  are  invariant 
with  respect  to  the  background  drift  rates  but  do  depend  on  the 
sampling  interval.  Note  that  for  a  fixed-difference  filter,  the 
number  of  “positive"  and  “negative"  samples  associated  with  the 
magnitude  response  is  one  higher  than  the  order  of  the  filter  (2  for 
a  first  difference  filter  in  Fig.  10.  for  example).  All  of  the 
preceding  statements  can  be  understood  by  examining  the  value  of 
the  coefficients  in  the  “tree  diagram"  of  Fig.  9. 


Target  in  View 
1 -Second  Sampling 


Some  of  the  limiting  assumptions  used  for  this  study: 

•  A  stochastic  quasi-linear  model  of  the  detector  process  is 
adequate. 

•  Noise  and  background  are  additive,  uncorrelated,  and 
generated  by  a  Gaussian  random  process  with  a  given 
power  spectral  density  (PSD). 

•  Target  is  additive. 

•  Sensor  and  filter  responses  are  linear. 

•  Background  and  target  velocities  are  constant  with  a 
fixed  detector  crossing  geometry. 

System  parameters  used  in  our  study  are  shown  in  Table  4. 


Fig.  7  shows  the  sampled  signal  magnitude  (sampling  interval 
of  1  sec)  and  noise  magnitude  spectra  up  to  a  frequency  range  of 
2  Hz  and  clearly  shows  the  effect  of  sampling  frequency  on 

Table  4.  Assumed  System  Parameters  for  Filter  Design 
to  Suppress  Background  Motion  Effects 
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Fig.  10.  First  Difference  Filter  Impulse  Response  (Same  for  all  Drift  Rates) 
l-Second  Sampling 

Fig.  1 1  shows  the  impulse  response  for  the  optimum  third- 
order  recursive  filter  which  was  derived  from  the  corresponding 
ideal  filter.  It  should  be  obvious  that  transversal  filter  approxi¬ 
mations  to  these  characteristics  would  require  many  more  than  just 
three  delay  stages.  The  impulse  responses  for  the  third-order 
optimum  recursive  filters  vary  with  both  the  sampling  interval  and 
the  background  drift  rate  since  the  filler  assumes  different 
coefficient  values  (given  in  Tables  2  and  depending  on  the 
sampling  interval  and  the  background  drift  rates.  Notice  also  that 
while  the  ideal  (non-realizable)  filter  is  non-causal,  the  optimum 
recursive  filter  being  realizable  has  a  causal  impulse  response.  Also, 
the  damping  of  the  impulse  response  depends  on  the  background 
motion.  For  the  optimum  recursive  filters  a  sharper  frequency 
domain  cutoff  and.  hence,  smaller  damping  coefficient  is  required 
for  high  background  drift  rate. 

Fig.  12  and  13  show  the  peak  normalized  magnitude  respon¬ 
ses  of  the  various  fixed-difference,  ideal,  and  third-order  optimum 
recursive  filters  to  a  frequency  range  of  1  Hz.  These  figures 
illustrate  the  advantages  of  optimum  recursive  over  fixed-difference 
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Fig.  12.  Magnitude  Responses  of  Various  Fixed  Difference  and  Third-Order 
Recursive  Filters  (liecond  Sampling  for  Optimum  Recursive  Filters) 
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Ftg.  14.  Effectiveness  of  Various  Filters  in  Suppressing  Background  Motion 
Effects  (l-Second  Sampling  Time) 
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Fig.  IS.  Effectiveness  of  Various  Filters  in  Suppressing  Background  Motion 
Effects  (2-Second  Sampling  Time) 


Fig.  13.  Magnitude  Responses  of  Third-Order  Optimum  Recursive  and  Fixed 
Difference  Filters  (2-Second  Sampling  for  Optimum  Recursive  Fillers) 

filters  of  first  through  fourth  orders.  For  a  sampling  time  interval 
of  I  sec.  lover  sampled  case)  Fig.  12  shows  that  for  high 
background  drift  rate  17.3  prad/sec),  the  third-order  recursive 
filter  has  a  greater  background  clutter  rejection  than  the  first-  and 
second-difference  filters  but  has  less  background  rejection  capabil¬ 
ity  than  the  fourth-difference  filter.  For  lower  background  drift 
rates  le.g..  1.2  prad/sec)  in  the  oversampled  case  (Fig.  12).  the 
third-order  recursive  filter,  optimized  for  Vj,  =  1 .2  prad/sec. offers 
a  broader  passband  than  the  second,  third-,  and  fourth-difference 
filters.  This  feature  allows  for  integration  against  the  white  shot 
noi.se  which  determines  the  SNR  over  the  background  drift  effects 
for  low  background  motion.  When  sampling  interval  coincides  with 
target  cell  dwell  time  as  in  Fig.  13.  the  advantage  of  recursive  filters 
over  fixed-difference  filters  is  again  observed  The  magnitude 
response  curves  of  the  two  third-order  recursive  filters  for  the 
extreme  values  of  the  background  drift  rates  span  the  magnitude 
response  curves  for  the  second-,  third-,  and  fourth-difference 
filten.  This  implies  that  by  adaptive  adjustment  of  its  coefficient 
values  in  response  to  the  existing  background  drift  rate,  the 
optimum  recursive  filter  can  indeed  perform  superior  to  the 
second-,  third-,  and  fourth-difference  filters  in  rejecting  back¬ 
ground  clutter. 

In  Fig.  14  and  15  are  plotted  the  maximum  peak  signal  to  rms 
noise  versus  the  background  drift  rates  of  interest  as  output  by  the 
various  fixed-difference,  ideal,  and  optimum  recursive  fillers. 


Figs.  14  and  15  assume  that  the  beginning  of  the  sampling  interval 
is  synchronized  with  the  passage  of  the  target  across  the  detector 
window,  such  that  during  the  sampling  interval  the  detector 
collects  the  maximum  possible  amount  of  the  target  energy. 

Fig.  14  and  15  show  the  effectiveness  of  the  optimum 
third-order  recursive  filters  (curves  labeled  TR  at  two  ends)  in 
suppressing  background  motion  effects  on  peak  SNR.  respectively. 
For  a  sampling  interval  of  I  sec.  Fig.  14  shows  that  for  lower  drift 
rates  (<4.9p  rad/sec),  the  third-order  recursive  filter  optimized  for 
the  same  drift  rates  have  more  than  a  factor  of  2  SNR  advantage 
over  equivalent  third-difference  filter;  at  higher  drift  rates  the  SNR 
improvement  is  le.ss  dramatic  as  the  effect  of  non-optimum 
sampling/integration  interval  affects  the  performance  of  the  recur¬ 
sive  filters  more  than  the  performance  of  third-  and  fourth- 
difference  filters.  Among  the  first-  through  the  fourth-order 
fixed-difference  filters,  for  all  ranges  of  drift  rates  and  sampling 
intervals,  the  performance  of  the  first-difference  filter  is  decidedly 
the  worst.  Sampling  time  interval  has  a  significant  influence  on  the 
relative  performance  of  various  orders  of  fixed-difference  filters. 
When  the  sampling  frequency  is  I  Hz  (Fig.  14).  for  lower  drift  rates 
«5.0  R  rad/sec),  the  second-difference  filter  is  superior  to  the 
third-difference  filter,  which  is  in  turn  superior  to  the  fourth- 
difference  filter.  For  the  range  of  drift  rates  considered  for  this 
study  (maximum  7.5  prad/sec).  there  is  no  advantage  in  using 
fourth-difference  filter  over  third-difference  filter  for  the  ove^ 
sampled  case,  but  as  shown  in  Fig.  1 5  and  Table  5.  for  higher  drift 
rates,  such  may  not  be  the  case.  Again,  for  the  range  of  low 
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Tables.  Summary  of  Resultt  Peak  Output  SNR  Gain  Calculations 
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backgroiiiKl  drill  considered,  lor  the  oversainpled  case,  the  peak 
SNR  performances  of  the  third- and  fourlh-dilTerence  filters  appear 
almost  unchanged  If'ig.  141.  as  the  shut  noise  is  unaffected  by  the 
drift  rate,  and  most  of  the  low  Vj,  clutter  has  been  eliminated  by 
high  difference  filters.  The  cross-over  of  the  third-order  recursive 
filler  mavimum  SNR  with  that  of  the  ideal  filler  in  Fig.  14  is 
caused  by  computational  appro.ximation  In  the  computer  simula¬ 
tion.  From  a  physical  standpoint  the  performance  of  all  filters, 
including  the  optimum  recursive,  suffers  with  increasing  back¬ 
ground  drill  rales. 

When  the  negative  effects  of  oversampling  are  eliminated  by 
using  a  sampling,  integration  Interval  equal  to  the  cell  dwell  time, 
the  superiority  of  the  optimum  recursive  filter  over  the  various 
orders  of  fixed-difference  filters  is  more  clearly  demonstrated.  This 
is  because  a  recursive  filter  can  compensate  for  improper  target 
integration  while  a  fixed-difference  filter  can  not.  While  a 
second-difference  filter  may  yield  higher  peak  SNR’s  (in  Fig.  15) 
than  third-  and  fourth-difference  filters  for  low  (<2.5  rad/sec) 
drift  rales,  for  higher  (>2.5  grad/see)  drift  rates,  the  higher  the 
difference  filler  order,  the  better  is  the  SNR  performance.  This  is 
because  lower  order  difference  fillers  degrade  the  target  less  for 
low  background  drift  rates.  For  higher  drift  rates  (>  7.0  g  rad/sec), 
the  third-order  recursive  filters  have  almost  a  factor  of  2  SNR  over 
the  third-difference  filter.  At  lower  drift  rales,  the  improvement  is 
less  noticeable  because  shut  noise  dominates  the  total  noise  input. 
Note  that  the  optimum  recursive  filter  SNR  performance  is  very 
close  I  within  I  O'?)  to  that  of  the  theoretically  "ideal"  filter  for 
drift  rates  shown  in  Fig.  15.  In  Fig.  15.  the  first-difference  filter 
has  peak  SNR's  of  less  than  I  for  drift  rates  higher  than 
.5.1  g  rad/sec.  It  is  not  surprising  that  recent  systems  have  not  used 
first-difference  filtering. 

Since  it  may  not  be  practical  to  adjust  the  filter  coefficients 
adaptively  for  each  drift  rate  for  an  optimum  recursive  filter,  it  is 
of  interest  to  study  the  behavior  of  a  recursive  filter  optimized  for 
a  particular  drift  rate  {usually  the  highest  expected  drift  rate)  for 
other  drift  rate  conditions.  Such  behavior  is  highlighted  in  Table  5 
for  a  third-order  recursive  filler  with  constant  filter  coefficients 
optimized  for  a  single  drift  rate  of  7, .5  grad/sec. 

Comparing  the  last  two  columns  of  Table  5.  we  see  that  much 
of  the  filler  performance  is  sacrificed  at  lower  drift  rates.  However, 
the  latter  filter  still  gives  higher  peak  SNR  over  the  range  of  drift 
rates  than  the  first-  through  the  fourth-difference  filters.  When  the 
effect  of  non-optimum  (I -sec)  sampling/integration  time  is  elimi¬ 
nated  from  the  problem  conditions,  we  again  find  by  comparing  the 
last  two  columns  of  Table  5  that  for  lower  drift  rates  (<3. 5  g  rad/ 
sec),  the  filter  optimized  for  7  .5  grad/sec  has  a  worse  SNR  than 
the  fixetWifference  fillers  except  the  first-order.  However,  as  the 
drift  rales  approach  the  drift  rate  for  which  the  filler  was 
optimized,  it  performs  better  than  any  of  the  fixed-difference 
filters. 


The  results  presented  here  show  that:  ( 1 1  the  use  of  optimized 
recursive  filters  instead  of  fixed-differencing  filters  can  result  in 
significant  improvements  in  SNR  performance  especially  at  higher 
(>5.0  g  rad/sec)  drift  rates;  (2)  in  addition,  by  oversampling  in 
time,  dropouts  due  to  phase  mismatch  between  the  target  and 
sample  time  can  be  lessened;  (3)  for  the  first-,  second-,  and 
third-difference  filters.  SNR  performance  is  slightly  improved  by 
oversampling  (sampling  interval  1/2  cell  dwell  time)  at  the  highest 
drift  rates  (6.0  and  7.3  g  rad/sec)  but  significant  SNR  degradation 
results  at  the  lower  drift  rates;  (4)  overall  performance  is  degraded 
by  not  using  an  optimum  target  signal  integration  time  and 
sampling  interval,  but  that  a  recursive  filler  can  be  adjusted  to 
compensate  to  a  large  extent  (filter  coefficients  can  be  set  in  a 
third-order  recursive  filter  so  that  both  differencing  at  low 
frequencies  to  attenuate  structured  background  and  integration  at 
high  frequency  to  attenuate  random  system  and  detector  shot 
noise  can  be  provided);  and  (5)  of  the  six  different  realizable  filters 
considered  in  the  study,  the  best  SNR  performance  was  achieved 
by  using  a  third-order  optimum  recursive  digital  filter  and 
oversampling  for  high  background  drift  rates  and  sampling  at  the 
target  cell  dwell  time  for  lower  drift  rates. 

IV.  Spatial  Filtering  Considerations 

Another  electro-optical  system  parameter  that  we  may  con¬ 
sider  varying  so  as  to  improve  the  SNR  performance  is  the  optical 
blur  circle  size.  To  gain  insight  into  the  performance  consequences 
associated  with  this  parameter,  we  have  posed  and  investigated  the 
following  problem: 

For  the  detection  of  moving  targets  immersed  in  a  drifting 
structured  background  with  a  staring  mosaic  type  sensor, 
derive  an  ideal  cascaded  pair  of  linear  spatial  and  temporal 
filters. 

For  this  problem  we  shall  assume  the  system  model  given  in 
Fig.  16.  In  addition,  we  shall  assume  the  following: 

1 .  The  total  input  noise  spectrum  can  be  written  as 

NtuT^-iat*  “  N[,(w  x)  +  N((u)j)  +  Njlujj) 

=  NbIcJx)  +  Njluij) 

where  Nlwx-w))  is  the  spectrum  of  the  total  input  noise. 
N^lwxl  's  Ih*-'  spectrum  of  the  background  spatial 
structure.  Nt(uit)  is  the  spectrum  of  the  background 
temporal  fluctualion.  and  Ns(wt)  is  the  spectrum  of 
temporal  random  noise  that  is  generated  in  the  detector, 
and  N,](w|)  is  the  total  temporal  noise. 

2.  The  phase  content  of  spectra  can  be  ignored  except  for 
the  computation  of  cross  spectral  densities. 
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Fig.  16.  Model  of  the  Caxatie  Spatial-Temporal  Filter 

3.  TIk-  backgrounil  is  drifting  with  velocity  V(,  and  the 
target  is  drifting  with  velocity  such  that  for  the 
background.  N(  =  ^b  •  “x  :  and  for  the  target  or  signal, 
"t  =  ^t  '  '^x 

4.  There  is  zero  correlation  between  the  inputs 

Nbtcjxl- iiftl 

The  ideal  temporal  filter,  which  will  maximize  at  its  output 
the  ratio  of  the  maximum  of  the  signal  to  the  rms  of  the  noise,  or, 
in  other  words,  the  SNR,  can  be  written  as* 

_  S(V,  ■  Jx>Hx(V,  •  Jx) 

HiltJ)  — ; - =i - =; - =; - ^ - - - -r-  (6) 

iNb'Vb  ■  “x>  Hx(Vb  ■  <Jx>l  -+  |Nd(‘et>l' 

Such  a  filter  has  the  property  that  the  greater  the  integrated  output 
spectrum  4,  ^ <u),)  d«t,  the  greater  will  be  the  output  SNR.  We 
can  thus  define  an  optimization  parameter  P  which  we  will  want  to 
minimize  as 


c(Z:7r‘‘“' 


Because  the  cross-spectra  terms  vanish.  Equation  7  can  be  reduced 
to* 

J  /lNb*Vb-^x>|-  |Hx(Vb-wxt|  - 
Jx'|-  iHxiVfWxll- 

d.,  (8) 

IStVfJx'l-  |Hx(V,  o;x)|  -  J 


From  this  equation  we  observe  that,  compared  to  no  spatial 
filtering  (Hx(urx)=  1 ).  a  spatial  filter  will  improve  the  signal  to 
baekground  stnicture  by  a  factor  of 


I  ^^x*^t  ’  ‘*’x*  I  ~ 

I  Hx(Vb  wxll  - 


and  the  signal  to  purely  temporal  noise  by  a  factor  of  |Hx(V(  U)x*|- 
at  each  frequency. 

Let  Gxlwj)  be  the  gain  in  signal  to  background  that  will  be 
achieved  by  a  spatial  lllter 


Gxlwf)  = 


_  |Hx(Vt  •  oJx)  I  ~ 


To  gain  insight  into  the  behavior  of  Gxlu>()  let  us  consider  the 
following  example.  Let  Nd(w{)  =  0  and  let  Hx(  w*x>  I’®  I*"-’  lowpass 
filter  of  Fig.  17.  Now  let  Hxtwx*  •I’'-’  highpass  filter  of  Fig.  18. 
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Fig.  17.  Illustration  of  the  Sigmi'tO'Background  Enhancement  Obtainabie 
From  a  Low-Pass  Spatial  Fitter 
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Fig.  18.  //tgA-Pau/Ttfrer 

It  appears  that  the  spatial  filter  will,  as  a  lowpass  filter,  yield  a 
gain  in  signal  to  background  as  a  function  of  frequency  for  a  target 
faster  than  the  background  drift.  A  highpass  spatial  filter  will  tend 
to  degrade  this  performance. 

Furthermore,  as  we  decrease  the  resolution  of  a  typical 
optical  system  we  should  observe  the  following; 

1 .  We  can  achieve  a  performance  gain  against  background  as 
expressed  in  Equation  9.  We  also  observe  as  shown  in 
Fig.  19  that  the  greatest  gain  will  occur  at  the  higher 
frequencies. 

2.  We  lose  in  performance  against  temporal  noise  by  a 
factor  G|( 

Gt(u) j)  =  Hx(V|  ■  Ux)  (10) 
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Fig.  19.  A  Typical  Optical  Low-Pass  Filter 

From  our  work  we  have  observed  that  as  we  increase  the  blur 
circle  size  we  always  lose  SNR  at  the  lower  drift  rates  but  we  gain  a 
potential  for  improved  SNR  at  higher  drift  rates.  Another 
observation  is  that  as  the  blur  circle  is  increased  the  breakeven 
point  for  SNR  versus  blur  circle  size  migrates  to  higlier  drift  rates 
so  that  for  each  drift  rate  there  will  be  some  single  optimum  blur 
circle  size  depending  on  the  drift  rate.  Another  observation  is  that 
as  the  blur  circle  is  increased,  temporal  filters  with  sharper  band 
edge  cutoffs  are  required  to  achieve  ideal  performance. 


I  o  achieve  such  sliarper  band  edge  cutoff,  as  shown  in  the 
previous  section,  one  must  use  more  generalized  recursive  filter 
structures  instead  of  just  simple  structures  associated  with  dif¬ 
ferencing  filters. 

V.  Summary /Conclusions 

Summarizing,  the  most  critical  problem  in  space-based  staring- 
mosaic -sensor  surveillance  systems  is  that  ideal  staring  operation 
cannot  be  achieved  with  non-synchronous  observer  orbits  looking 
at  earth  backgrounds,  background  appears  to  drift  everywhere  in 
the  I  0\'  eveept  the  stabilized  center.  The  objectives  of  this  study 
have  been  threefold: 

1 .  1  o  characterize  earth  background  drift  versus  orbit. 

2.  Compare  effectiveness  of  signal  processing  algorithms  in 
rejecting  the  background  clutter  noise. 

-V  I'se  of  spatial  filtering  effects  in  suppressing  background 
clutter. 

the  following  were  the  most  significant  results  of  our  study; 

1.  I  or  lived  angular  FOV: 

a.  Background  drift  in  angular  measure  (arc-sec/sec) 
decreases  markedly  with  observer  altitude  (Fig.  4a). 

b.  Angular  drift  rate  varies  considerably  with  viewing 
geometry.  A  factor  of  2.5  to  1.5  increase  results 
(depending  on  orbit  altitude)  from  increasing  LOS 
zenith  angle  (angle  between  LOS  and  local  vertical 
at  LOS  intersection  with  the  ground)  from  30  to 
80°  (Fig.  4b).  Since  range  also  increases  with 
increasing  zenith  angle,  apparent  linear  drift  on  the 
ground  increases  even  more. 

c.  Background  drift  in  linear  measure  (meters/sec) 
changes  little  with  altitude  (Fig.  6). 

2.  Assuming  drift  varies  linearly  with  distance  from  FOV 
center,  average  angular  drift  in  a  circular  FOV  is 
two-thirds  of  the  worst-case  drift  at  the  edge  of  the 
FOV.  Average  drift  at  the  periphery  of  a  circular  FOV  is 
approximately  1Q7<  of  the  maximum  drift  at  one  or 
more  corners  of  the  circumscribed  square  FOV.  For 
fixed  coverage  on  the  ground,  requiring  smaller  FOV  for 
higher  observer  altitude,  linear  drift  (m/sec)  and  hence, 
background  clutter  noise  decreases  with  altitude:  this 
implies  higher  orbits  should  be  favored  for  aircraft 
detection. 

3.  One  can  obtain  significantly  better  SNR  performance 
from  an  optimized  recursive  filter  than  from  an 
equivalent  fixed-differencing  filter,  especially  at  higher 
drift  rates.  For  equal  performance.  FOV  area  for  a 
third-order  recursive  filter  equals  four  times  the  FOV 
area  for  a  third-difference  filter:  the  latter  implying 
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step-stare  update  rate  of  four  times  as  much  as  for  a 
recursive  filter  (Fig.  14  and  15).  Higher  order  filters  are 
quite  effective  in  improving  clutter  noise  rejection.  In 
one  sample  case,  for  higher  drift  rates,  a  third-order 
recursive  filter  permitted  only  one-half  of  the  clutter 
leakage  of  a  third-differencing  filter  which  in  turn 
permits  about  one-half  of  the  clutter  leakage  of  a 
second-difference  filter. 

4.  While  the  performance  of  fixed-differencing  filters  is 
severely  degraded  (if  the  sample  rate  is  higher  than  the 
target  cell  dwell  time),  recursive  filters  can  be  adjusted  to 
compensate  for  the  non-optimum  detector  integration 
time.  In  addition,  by  oversampling  in  time,  dropouts  due 
to  phase  mismatch  between  the  target  and  the  sample 
time  can  be  minimized.  Parameters  can  be  set  in  a 
third-order  recursive  filter  so  that  both  differencing  at 
low  frequencies  to  attenuate  structured  background  and 
integration  at  high  frequency  to  attenuate  random 
system  and  detector  shot  noise  can  be  provided.  To 
summarize,  we  have  shown  that  superior  performance 
can  be  achieved  using  an  optimum  recursive  digital  filter 
(and  oversampling  in  time  for  high  background  drift 
rates  and  sampling  at  target  cell  dwell  time  for  low  drift 
rates)  than  from  a  fixed  differencing  filter  operated  at 
one  sample  per  target  cell  dwell  time. 

5.  A  recursive  filter  that  is  optimum  for  a  given  background 
drift  rate  will  exhibit  somewhat  higher  SNR's  than  fixed 
difference  filters  when  the  background  drift  rate  is 
actually  lower  than  assumed  (Table  5). 

6.  Peak  signal  to  rms  background  clutter  can  be  improved 
by  increasing  the  optical  blur  circle  size  and  using  a 
temporal  filter  with  a  sharper  band  edge  cutoff.  How¬ 
ever.  as  the  blur  circle  is  increased,  peak  signal  to  rms 
shot  and  system  noise  will  be  decreased. 
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